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In the present work, the heat transfer study focuses on assessment of the impact of bed temperature on the local 


heat transfer characteristic between a fluidized bed and vertical rifled tubes (38mm-O.D.) in a commercial circu- 
lating fluidized bed (CFB) boiler. Heat transfer behavior in a 1296t/h supercritical CFB furnace has been analyzed 
for Geldart B particle with Sauter mean diameter of 0.219 and 0.246mm. The heat transfer experiments were 


conducted for the active heat transfer surface in the form of membrane tube with a longitudinal fin at the tube 


crest under the normal operating conditions of CFB boiler. A heat transfer analysis of CFB boiler with detailed 


consideration of the bed-to-wall heat transfer coefficient and the contribution of heat transfer mechanisms inside 


furnace chamber were investigated using mechanistic heat transfer model based on cluster renewal approach. The 


predicted values of heat transfer coefficient are compared with empirical correlation for CFB units in large-scale. 
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Introduction 


The issue of heat transfer process between the flui- 
dized bed and the boiler surface in the form of water 
membrane wall plays a critical role in the case of circu- 
lating fluidized bed reactors used for the heat recovery 
and the clean or high efficiency combustion of fossil fu- 
els. Particular attention has been recently drawn to the 
reduction of the energy intensity of the energy generation 
process, while meeting the environmental protection re- 
quirements laid down by the IED Directive of the Euro- 
pean Parliament and the Council [1]. Understanding the 
complex gas-solid flow pattern and heat transfer pro- 
cesses occurring during solid fuel combustion under cir- 
culating fluidized-bed conditions constitutes a factor de- 
cisive to the quality of the final product obtained. This 


the heat exchangers, while taking into account the loca- 
tion and geometry of the heatable surfaces. Understand- 
ing the mechanisms of heat flow, being in a close corre- 
lation with the prevailing hydrodynamic conditions, will 
provide the necessary knowledge of the process scaling 
and the design of heat exchange surfaces. The correct 
sizing of active heat transfer surfaces inside CFB furnace 
ensures proper operation, load turndown and also affect 
furnace temperature. The furnace temperature is of im- 
portance for the formation and reduction of harmful mat- 
ter, which influences the economic indices of commercial 
CFB boilers. 

The aims of the current work were to: (i) experimental 
investigation of bed temperature and suspension density 
profiles in a 1296t/h supercritical CFB boiler; (ii) present 


applies to the method of receiving and supplying heat to 


Nomenclature 
a coefficient in Eq. (2), - 
b coefficient in Eq. (2), - 


CFB circulating fluidized bed, - 


Ung minimum fluidization velocity, m/s 
U, terminal velocity of bed particles, m/s 
Q boiler load, % 


c coefficient in Eq. (2), - 

Cg specific heat of gas, kJ/(kg:K) 

Cp specific heat of bed particle, kJ/(kg:K) 

D; hydraulic diameter, m 

dp mean bed particle size, mm 

e emissivity, - 

f fractional of wall covered by clusters, - 

Gsh lateral solid circulation flux, kg/(m’s) 

g acceleration due to gravity, m/s? 

H furnace height, m 

h bed-to-wall heat transfer coefficient, W/(m°K) 

he cluster heat transfer coefficient, W/(m°K) 

hg gas convection heat transfer coefficient, W/(m’K) 

n particle convection heat transfer coefficient, 
d W/(m'K) 

hy, wall contact heat transfer coefficient, W/(m’K) 

k thermal conductivity, W/(m:K) 

Le cluster characteristic travel length, m 

MCR maximun continuous rating, - 

p pressure, kPa 

Ap pressure drop, kPa 

Pr Prandtl number, - 

T temperature, K 

te cluster residence time, s 

U. cluster descent velocity, m/s 

U; superficial gas velocity, m/s 


overall heat transfer coefficient distributions; and (iii) 
depict the contribution of heat transfer mechanisms in- 
side furnace as a function of bed temperature. 

In summary, this work addresses an existing gap in the 
literature data regarding the modelling of and experi- 
mentation of heat transfer to tubes in furnace chamber 
[2-6], by providing experimental data from CFB unit in a 
large-scale. The high labour intensity and cost of carrying 
out large commercial-scale tests are the cause of the lack 
of publications of the results of investigations of heat 
transfer processes on a large commercial scale [7-10]. In 
this situation, carrying out comprehensive studies would 
make it possible to establish the optimal heat transfer 
conditions in correlation with the bed hydrodynamics. 


Heat Transfer Process In CFB Furnace 


Heat transfer between the fluidized bed and containing 
walls involves three modes: (i) convection by particles, 
(ii) convection by gas, (iii) radiation from clusters, and 
also (iv) radiation from dispersed phase. Contribution of 
individual heat transfer mechanisms are considered sep- 


fraction of particles in the dispersed 


: phase, - 
Z height above air distributor, m 
Az distance between pressure taps, m 
Greek letters 
ô wall layer thickness, m 
Of fin depth, mm 
`z gas gap thickness, m 
E cross-sectional bed average voidage, - 
u viscosity, kg/(m:s) 
p density, kg/m’? 
Pb suspension density, kg/m? 
(o Stefan Bolzmann’s constant, W/(m K$) 
Superscripts 
ad air dried basis 
ar as received 
Subscripts 
b bed 
c cluster phase 
cr radiation from cluster phase 
dr radiation from dispersed phase 
g gas 
max maximum value 
p particle 
w wall 
00 core 


arately due to the complexity of the bed hydrodynamics. 
Mechanistic models based on cluster renewal approach 
commonly assume that the walls are not covered by a 
continuous stream of solids. In the case of CFB boilers 
each portion of the active heat transfer surface (i.e. water 
membrane wall) is intermittently washed by down- 
flowing clusters of bed particles and dispersed medium. 
Thus, one part of the heat transfer surface is in contact 
with cluster while the rest is in contact with the gas or 
dispersed phase. This fact is due to a major feature of 
core-annulus flow structure inside CFB furnace as indi- 
cated by Fig. 1. 

So, the bed-to-wall heat transfer coefficient can be 
stated in the following equation (1) as the sum of heat 
transfer mechanisms contributions [12-14]: 

h= fhy+(1-f\hg + fhet- fas A) 

where the subscripts p, g, rc and rd denote particle 
convection, gas convection, radiation from clusters and 
radiation from dispersed phase, respectively. In equation 
(1), the parameter f represents the fractional wall cover- 
age by clusters and is estimated using the hydraulic di- 
ameter and height of the CFB furnace in a non-dimen- 


sional form: 


f=1-expf-a(l—e)"(D,/#1)'| 2) 
To derive values of the coefficient a, b, c in equation 
(2) heat transfer data from four CFB boilers in large-scale 
with a different equivalent diameter 1.6 m (12MW»n), 5.2 
m (20MW,), 6.2 m (145MW,,), 10.6 m (170MW.) and 
also at different elevations above the secondary air 
supply: 11.5 m, 25.0 m, 26.0 m and 30 m were ana- 
lyzed and then the coefficients have been estimated using 
multiple regression analysis [9]. 
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Fig. 1 Single cluster formation and gas gap in the vicinity of 
the water membrane wall inside CFB furnace [11]. 


In the cluster renewal model of heat transfer for CFB 
boilers, particle convection heat transfer coefficient com- 
prises unsteady heat conduction into clusters A. and the 
conduction across the thin gas layer between the cluster 
and the heat transfer surface h,,. Finally, the particle con- 
vection heat transfer coefficient h, is given by: 


-1 -1 
d,ô 
(EE) - Te i). @) 
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Formulas needed to estimate the physical and thermal 
properties for the clusters and gas have been detailed 
described by Blaszczuk et al. [12] and will not be re- 
peated in this work. 

When the vertical heat transfer surfaces (i.e. mem- 
brane walls) exposes to fluidized bed at large Archimedes 
number or dilute lean phase, the gas convection compo- 
nent in overall heat transfer coefficient may be expressed 
as follows [15]: 


7 037, \021 
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where k, denotes thermal conductivity of the gas, d, 
represents mean bed particle size, U, means terminal ve- 


locity of the bed particles, c, and c, are specific heat of 
gas phase and particle, respectively. In the above equa- 
tion (4), g is acceleration due to gravity, p, represents 
particle density and p4 denotes dispersed phase density. 
In the current heat transfer study, dispersed phase density 
is found by the following equation: 

Pa = PpY + pg(l-Y). (5) 

In this correlation, Y is the volumetric concentration of 
bed particles in the dispersed phase and the value of this 
parameter is recommended as 0.001% [16]. In the heat 
transfer model, the dispersed phase contains a small 
concentration of bed particles (i.e.dilute gas-particle 
mixture) and gas. 

The radiation heat transfer from the cluster to the wa- 
ter membrane wall is similar to that between two infinite 
parallel plates with a view factor 1. Thus, the cluster rad- 
iation component of the bed-to-wall heat transfer coeffi- 
cient is calculated from the expression: 


o(ré ý) 
RS (T-T, We segla) 


where o is the Stefan-Boltzmann constant, T„ repre- 
sents active heat transfer surface temperature, T, means 
cluster temperature which is estimated from the empirical 
correlation for radial temperature distribution proposed 
by Golriz [17]. In above equation (6), ec and e, are the 
emissivities of the cluster and wall, respectively. The 
membrane wall emissivity e,, is assumed to have a con- 
stant value 0.8, as for oxidized steel. For analysis of heat 
transfer in commercial CFB combustors, the cluster 
emissivity can be derived by equation (7) which takes 
account of bed particle radiation: 

e, =0.5(I+e,). (7) 

Here e, denotes particle emissivity and is equal to 0.7 
as suggested Basu [18]. Thus, the cluster emissivity is 
larger than particle emissivity. 

At high furnace temperature (>973 K) and also region 
with low suspension density (<5 kg/m?) inside furnace 
chamber of CFB boilers [18-20], the dispersed radiation 
component of the overall heat transfer coefficient be- 
comes important and is estimated using the following 
equation: 


(6) 


TE = (8) 
(7, -Tez +e," )-1) 
where T, denotes bed temperature, eg is the emissivity 
of dispersed phase which may be calculated by the 
Brewster’s relationship for isotropic scattering [21]. The 
expression for e, is given as: 


e 
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For isotropic scattering, B can be taken to equal about 
0.5. Other denotations in equation (8) are explaned in 


equation (6). 
Performance Tests and Experimental Conditions 


The heat transfer experiments were performed in a 
1296 t/h supercritical CFB boiler located in Tauron Gen- 
eration S.A. Lagisza Power Plant, Poland. More informa- 
tion about the CFB facility can be found in the works [12, 
13, 22, 23], in which particularly is presented construc- 
tion data, cross-sectional area, arrangement of measuring 
points and also data acquisition system. 

The planning of large commercial-scale experiments 
required a series of preliminary tests to be performed, 
which involving the observations of the variability of 
process parameter change rate. Based on these observa- 
tions, the frequency of process data archiving in accor- 
dance with the Shannon-Kotelnikov theorem [24] and the 
duration of an individual test were selected. 

To assessment of the impact of bed temperature on the 
local heat transfer characteristic between a fluidized bed 
and vertical rifled tubes, experiments were conducted 
during a four-day period. 

Each test lasted eight hours at four different boiler 
loads (i.e. 40%MCR, 60%MCR, 80%MCR and 
100%MCR). The eight hours for each test were needed in 
order to establish the repeatability of experimental results. 
The measurement of furnace data (i.e. static pressure and 
bed/wall temperature data) was carried out under stable 
operating conditions. During performance tests on the 
supercritical CFB boiler with capcity 966MW, the mass 
flow of feeding materials supplied into furnace chamber 
(i.e. coal and limestone) was kept on constant level. 

Air fossil fuel firing conditions of performance tests 
are summarized in Table 1. 

For the performance tests described in this paper, bed 
temperatures between 1063 K and 1183 K were achieved 
by burning Polish bituminous coal (Ziemowit coal mine, 
Poland) with approximately 20% excess air. The characte- 
ristic of fossil fuel used in the performance tests are given 
in Table 2. 

Determination of proximate parameters of the bitu- 
minous coal was made in accordance with applicable 
standards for fossil fuel in Poland. Meanwhile, the ulti- 
mate analysis data were obtained by means of the LECO 
TrueSpec™ analyzer. All analytical data for fossil fuel 
are averages from five repetitions for each fuel constituent. 
Table 1 Operation parameters of performance tests 


Parameter Unit Overall range 
Boiler loads, Q %MCR 40-100 
Superficial gas velocity, U, m/s 3.89-5.11 
Minimum fluidization velocity, Umf m/s 0.0164-0.0243 


Voidage at minimum fluidization 
. - 0.40-041 
velocity, Enf 


Solids circulation rate, Gs kg/(m’s) 21.4-25.5 
Sauter mean particle diameter, d32 mm 0.219-0.246 
Suspension density, p, kg/m? 1.55-6.14 
Bed temperature, Tp K 1063-1183 
Wall temperature, Tw K 636-735 
Pressure drop, Ap kPa 6.86-8.25 


Table2 Proporties of bituminous coal used in the tests 


Ultimate analysis 


(it dted basis) Unit Overall range 
C™, carbon wt.% 52.32-56.50 
H“, hydrogen wt.% 4.02-4.74 
oO, oxygen wt.% 5.35-6.98 
N“) nitrogen wt.% 0.73-0.84 
S% sulphur wt.% 0.87-1.21 
Proximate analysis (as-received) 
Q", Caloric value MJ/kg 20.11-22.79 
V“, Volatile matter wt.% 25.61-30.37 
A™, Ash wt.% 9.10-20.11 
M”, moisture wt.% 11.81-19.13 


In supercritical CFB boiler above the secondary air in- 
jection level (transport zone) are located the active heat 
transfer surfaces i.e. membrane wall, wing wall evaporators 
and also radiant super heater. In the current study is consi- 
dered heat transfer at membrane wall inside furnace cham- 
ber. The membrane wall is the main absorbing heat transfer 
surface from products of solid fuel combustion process. 
The membrane walls are constructed in the form the mem- 
brane rifled tubes with a longitudinal fin at the tube crest. 
The veritical rifled tubes have a 38mm outside diameter 
with spacing of 63mm. More detailed information on 
geometry of membrane walls can be found in the follow- 
ing publications [14, 25]. 


Results and Discussion 


Experimental data used in this heat transfer study are 
given in dimensionless scale by reason of confidential in- 
formations about industrial CFB boilers in large-scale. Be- 
sides, suspension density was normalized by maximum 
value obtained during all tests. 

In this heat transfer study, the concept of mechanistic 
heat transfer model is used to investigate the heat transfer 
mechanism behavior near the vertical rifled tubes, ac- 
cording to the cluster renewal approach. The bed-to-wall 
heat transfer coefficient and contribution of heat transfer 
components along furnace chamber are generated based on 
the experimental conditions given in Table 1. 

Bed-to-wall heat transfer coefficient distribution along 
furnace height at different bed temperature for a 1296 t/h 


supercritical CFB boiler is depicted in Fig. 2. 
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Fig. 2 Bed-to-wall heat transfer coefficient variation with fur- 
nace height of CFB boiler at different bed temperature. 


With increase in bed temperature, heat transfer coeffi- 
cient increases for the same non-dimensional distance Z/H 
from the air dirstributor level. The bed-to-wall heat transfer 
coefficient increases with bed temperature for higher frac- 
tional of wall covered by clusters. This fact is due to the 
combined effect of increased cluster heat transfer and 
reduced gas gap resistance due to higher gas thermal 
conductivity. Other reason for the higher bed-to-wall heat 
transfer coefficient arises from suspension density varia- 
tions along furnace height of CFB boiler. 

Figure 3 illustrates a strong influence of suspension 
density on the bed-to-wall heat transfer coefficient. 

In our heat transfer study, solid suspension density 
was calculated by the following empirical correlation: 

1 Ap 

g Az” 

where Ap represent pressure drop and Az denotes dis- 
tance between pressure taps. Attrition and gas-particle 
acceleration effects were neglected in equation (10). 

The bed temperature has an essential effect on the heat 
transfer data at high relative suspension density. In the 
case of low normalized suspension density, the bed tem- 
perature has negligible impact of the overall heat transfer 
coefficient. In addition, it can be seen from Fig. 3 that the 
influence of solid suspension density diminishes with bed 
temperature. In Fig. 3, the heat transfer data obtained 
using the cluster renewal approach were correlated by the 
non-linear function using the regression analysis. The 
heat transfer grey curves in Fig. 3 have steeper slope in 
the experimental data region. In addition, it can be seen 
from Fig. 3 that there is not much difference between the 
experimental data points for the cluster renewal approach 
and the published experimental data for Chalmers and 
Dalhousie boiler [9, 26] at the suspension density varied 


p= (10) 


in the range of 2-4 kg/m’. In the case of suspension den- 


sity covering the range of 5.44-6.14 kg/m’ was observed 
significant difference. 
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Fig. 3 Bed-to-wall heat transfer coefficient distribution as a 
function of the bed temperature and the normalized sus- 
pension density. 


Figures 4-7 illustrate contribution of heat transfer me- 
chanisms at different bed temperature inside CFB furnace. 
Contribution of heat transfer mechanisms is drawn as indi- 
vidual colour and type line for each heat transfer mode. In 
this heat transfer study, the gas convection contribution was 
not significant in the bed-to-wall heat transfer coefficient, 
as demonstrated by Figs 4-7. 
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Fig. 4 Contribution of heat transfer mechanisms as a function of 
bed temperature at Z/H=0.25 


At non-dimensional height coordinate Z/H of 0.25 (low- 
er region of transport zone) was recognized high bed par- 
ticles concentration and more clusters in comparison with 
other elevations above the fluidization grid (i.e. 0.5<Z/H< 
0.87). Due to high suspension density in the vicinity of the 
vertical riffled tubes, particle convection component in heat 
transfer mechanisms prevails in the overall heat transfer 
coefficient. The high particle convection heat transfer re- 
sults from the thermal conductivity of cluster. The cluster 
thermal conductivity increases with bed temperature due to 
higher gas thermal conductivity. The resistance due to gas 
gap thickness between cluster and riffled tubes was reduced 
due to increased gas thermal conductivity. The value of 


h,a'h increases with the increase of bed temperature. The 
experimental tests show that radiation from cluster phase 
not depend upon furnace temperature direct above the slope 
section of furnace chamber. 

At the half CFB furnace height (Z/H=0.5), radiative heat 
transfer plays dominant role in heat transfer mechanisms 
whereas bed temperature exceed 1125 K. This fact indi- 
cates significant role of suspension density and also furnace 
temperature (i.e. wall and bed temperature) on contribution 
of heat transfer modes between the fluidized bed and the 
active heat transfer surface inside CFB furnace. The slope 
of particle convection curve is directly propotional to bed 
particle concentration and bed temperature. The same trend 
was also confirmed by [12]. The radiation from clusters 
slowly decreases with the increase of radiation from dis- 
persed phase. 
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Fig. 5 Contribution of heat transfer mechanisms as a function of 
bed temperature at Z/H=0.5 


The variations of heat transfer mechanisms contribution 
in the exit region of furnace chamber occupied by mem- 
brane wall and radiant superheater SH II (transport zone at 
0.65<Z/H<0.87) are presented in Figs. 6-7. 

The shape of profiles for A/h, h,./h and also h,q/h con- 
tribution is almost the same, both for Z/H=0.65 and Z/H= 
0.87. The upper part of furnace chamber was operated in a 
dilute phase (p,<1.55 kg/m*) whereas the middle and lower 
regions of transport zone were a dense phase (p;> 6 kg/m’). 
As can be seen in Figs. 6-7, dilute phase radiation heat 
transfer h,a increases with bed temperature and is characte- 
rized by a linear function for non-dimensional distance 
0.65<Z/H<0.87. Neverthelees, the radiation from dispersed 
phase at 0.25<Z/H<0.5 is characterized by exponential 
function. This conclusion is consistent with the results pre- 
sented in Figs. 4-5. Moreover, a similar trend of the heat 
transfer components variation was observed by Basu [18], 
when CFB unit was operated at 70% load condition. Due to 
sufficient low suspension density (dilute phase) and high 
bed temperature in the exit region of furnace chamber, rad- 
iation from dispersed phase has dominant role in heat 


transfer mechanisms from core region to vertical riffled 
tubes. 
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Fig. 6 Contribution of heat transfer mechanisms as a function of 
bed temperature at Z/H=0.65 
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Fig. 7 Contribution of heat transfer mechanisms as a function of 
bed temperature at Z/H=0.87 


The performance tests show that there is interplay of bed 
particle concentration and clusters in contributions of par- 
ticle convection and radiation from clusters. Under dilute 
phase conditions the active heat transfer surface was cov- 
ered by less clusters (ca. 13%) than in a dense region of 
transport zone (ca. /=68%). 


Conclusions 


In this heat transfer study, the mechanistic heat trans- 
fer model with integrates bed hydrodynamics in CFB 
furnace was used to assessment of the effect of bed tem- 
perature on heat transfer behavior in CFB coal combustor. 
A heat transfer analysis was carried out between fluidized 
bed and vertical riffled tubes based on cluster renewal 
approach. Heat transfer data were compared with pub- 
lished literature data for CFB units in large-scale. The 
experimental data agree quite well with the results for 
CFB boilers in large-scale. The experimental measure- 
ments of furnace data (bed temperature and suspension 


density) show that bed-to-wall heat transfer coefficient 
was strongly correlated with furnace height. The bed 
temperature and suspension density affected the overall 
heat transfer coefficient. This conclusion was also con- 
firmed by authors [18, 19, 28]. Obtained heat transfer 
data prove the effect of fraction of wall exposed to clus- 
ters on particle heat transfer coefficient, especially at 
non-dimensional distance of 0.25 and also 0.5 from flui- 
dization grid. The low particle concentration (suspension 
density) and low value of the parameter f contribute to 
decrease of convection heat transfer mechanism. Mean- 
while, the radiation from dispersed phase proportion in- 
creased with increase in bed temperature. In the upper 
part of combustion chamber occupied by superheater SH, 
the contribution of the radiative heat transfer coefficient 
was dominant in the bed-to-wall heat transfer. This is due 
to sufficient high bed temperature and also low suspen- 
sion density. The contribution of radiation from clusters 
does not exceed 20% and 12% for 0.25<Z/H<0.5 and 
Z/H>0.65 respectively. The gas convection heat transfer 
mechanism was negligible during all performance tests. 
The heat transfer data will be useful in the design and 
scale-up of heat exchangers operated in circulating flui- 
dized bed conditions. Very important and useful for the 
development of computation methods and the design of 
commercial CFB boilers is the comprehensive study of 
the impact of bed temperature on the heat transfer me- 
chanisms behavior. The observed variation of heat trans- 
fer mode will give fundamental design basic for fluidized 
bed systems 
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